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The exact momentum conservation laws for the nonlinear gyrokinetic Vlasov-Poisson equations 
are derived by applying the Noether method on the gyrokinetic variational principle [A. J. Brizard, 
Phys. Plasmas 7, 4816 (2000)]. From the gyrokinetic Noether canonical- momentum equation derived 
by the Noether method, the gyrokinetic parallel momentum equation and other gyrokinetic Vlasov- 
moment equations are obtained. In addition, an exact gyrokinetic toroidal angular-momentum 
conservation law is derived in axisymmetric tokamak geometry, where the transport of parallel- 
toroidal momentum is related to the radial gyrocenter polarization, which includes contributions 
from the guiding-center and gyrocenter transformations. 
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^ ' I. INTRODUCTION 

Pk 

The property of exact energy conservation has long been a stringent test on the accuracy of nonlinear gyrokinetic 
numerical simulation codes [lj . While the gyrokinetic energy conservation law was initially constructed directly from 
the nonlinear gyrokinetic equations (2j-|5[, a more general approach relies on the existence of a variational formulation 
for the nonlinear gyrokinetic equations frjl-floj and the use of the Noether method [ll[ associated with the symmetry 
^ ' of the gyrokinetic Lagrangian density with respect to infinitesimal time translations. 

The turbulent transport equations for parallel momentum and toroidal angular momentum in magnetized plasmas 
have recently been the focus of intense investigation These gyrokinetic momentum transport equations are 

generally derived as moment equations of the nonlinear gyrokinetic Vlasov equation [l3l - fT7l | . The present paper is 
. concerned with an alternative derivation of the momentum and angular momentum conservation laws for the nonlinear 
Oh' gyrokinetic Vlasov-Poisson equations by the Noether method, which associates a conservation law to each symmetry 
of the gyrokinetic Lagrangian density with respect to infinitesimal space-time translations and rotations. The energy- 
' momentum conservation laws associated with the guiding-center Vlasov-Maxwell equations have previously been 
investigated with a variational formulation by Pfirsch and Morrison [ljj]. The present work generalizes Pfirsch and 
l/") ■ Morrison's work [l8| to the gyrokinetic Vlasov-Poisson equations, first, and the gyrokinetic Vlasov-Maxwell equations 
in a future paper. 

The gyrokinetic Noether canonical-momentum equation derived for the gyrokinetic Vlasov-Poisson equations is 
an exact equation that describes the transport of gyrokinetic canonical momentum associated with the gyrocenter 
particles (represented by the gyrocenter Vlasov distribution) and the electrostatic fluctuations (represented by the 
gyrokinetic Poisson equation), in the presence of a nonuniform background magnetic field. While momentum can be 
exchanged between the gyrokinetic Vlasov-Poisson fields and the background magnetic field when the latter does not 
possess exact spatial symmetries (e.g., toroidal and poloidal directions in stellarator geometry (i~9l - [2l1 |). the canonical 
toroidal angular momentum carried by the gyrokinetic Vlasov-Poisson fields is conserved exactly in an axisymmetric 
tokamak magnetic field, which allows the transfer of momentum between the gyrocenter Vlasov distribution and the 
electrostatic-field fluctuations. From the gyrokinetic Noether canonical-momentum equation derived by the Noether 
method, the gyrokinetic parallel momentum equation and other gyrokinetic Vlasov-moment equations are obtained. 

The remainder of the paper is organized as follows. In Sec.|lU we review the variational derivations [til. l9l. [Io| of the 
nonlinear gyrocenter canonical and noncanonical Hamilton equations and the nonlinear gyrokinetic Vlasov-Poisson 
equations. In Sec. lIIIl we apply the Noether method (briefly reviewed in App.lA"]) to derive the exact gyrokinetic linear 
and angular momentum conservation laws. In Sec. lIVl we present an application of the gyrokinetic angular momentum 
conservation law by deriving the parallel-toroidal gyrokinetic momentum conservation law in axisymmetric tokamak 
geometry. Lastly, we summarize our work and discuss future work in Sec. fVl 



II. GYROKINETIC VARIATIONAL PRINCIPLES 



The variational derivations of the nonlinear gyrokinetic Vlasov-Poisson and Vlasov-Maxwell equations and the 
gyrocenter Hamilton equations guarantee that these equations possess important conservation laws such as energy 
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and momentum. The purpose of this Section is to review these variational principles before they are used in the 
remainder of the paper. 

A. Gyrocenter Hamilton variational principle 

The single-particle variational principle used for the derivation of the gyrocenter Hamiltonian dynamics of a charged 
particle of mass m and charge e is based on the noncanonical phase-space Lagrangian [3, HJ 

T gy = (- A + p|| b) -dX + (mc/e) /idC - H &y dt = r « dz ° - H sy dt > i 1 ) 

where the background magnetic field B = V X A is defined in terms of the time- independent magnetic vector potential 
A, while the gyrocenter phase-space coordinates are the gyrocenter position X, the gyrocenter parallel momentum 
pu , the gyrocenter magnetic moment and the gyrocenter gyroangle £. The gyrocenter Hamiltonian H gy is defined 
as 

#gy( x >P||,jM; 0i ) = #gc(X,p||,/i) + e$gy(X,/x,t), (2) 

where H gc = /iB + p^/2m denotes the (unperturbed) guiding-center Hamiltonian. The gyrocenter potential $ gy is 
defined in terms of the first-order fluctuating guiding-center electrostatic potential 

0i gc = T~Vi = 0i(X + p gc ,f), (3) 



$ gy = e^ige) - y 0igc} g \ + •■■ , (4) 

where p gc denotes the guiding-center gyro radius (which contains higher-order contributions due to the nonuniformity 
of the background magnetic field I22I l23|Q (• • • } denotes the gyroangle- averaging operation, and the guiding-center 
Poisson bracket { , } gc is defined as |22l. l23j 

f „ flfdFdG dF dG\ B* f^^dG dF „\ cb „ 

where 

B* e VxA* = B + (cp||/e)Vxb = Sjf b + p\\ — X (b-Vb) (6) 

and -Bj| = b • B* = B + (cp\\/e) b • V X b. The second-order ponderomotive potential in Eq. §4$ is expressed in terms 
of the first-order function Si , which is formally defined as the solution of the first-order equation [2J, |6| 

^L-I = e 0i gc ee e ^ lgc - (0i gc )) , (7) 

where d gc /dt = d/dt + { , -ff gc } gc denotes the unperturbed (guiding-center) Vlasov operator. To lowest order in the 
gyrokinetic space- time-scale ordering [f|, the solution for S\ is approximated as 

Si s ^ / 0i gc dC, (8) 

where higher-order terms (not considered in standard gyrokinetic theory Q) have been derived in Ref. 24]. We note 
that, while the gyroangle-independent potential (4>i gc ) contributes to the linear (first-order) perturbed gyrocenter 
Hamiltonian dynamics, the gyroangle-dependent potential </>i gc contributes to the (second-order) gyrocenter pondero- 
motive Hamiltonian in Eq. @ as well as gyrocenter polarization effects [see Eq. (|5T)l below] . Lastly, high-order terms 
in Eq. (|4]) have recently been derived at third order in Ref. 25], which yield nonlinear gyrocenter polarization effects. 
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1. Gyrocenter noncanonical Euler- Lagrange equations 

The noncanonical gyrocenter equations of motion are obtained from the variational principle 

- */r„ . / M- ( Uo!J 4^ - ) «, (9) 

which, for arbitrary virtual displacements SZ a in gyrocenter phase space (that vanish at the integration boundaries), 
yields the noncanonical gyrocenter Euler-Lagrange equations 

- V = s#' (10) 

where the antisymmetric Lagrange matrix w Q( 3 = daT^ — c^Lq, is derived from the symplectic part of the phase-space 
Lagrangian ([lj, while the gyrocenter Vlasov operator d gy /dt = d/dt + { , H gy } gc is defined in terms of the gyrocenter 
Hamiltonian @ and the guiding-center Poisson bracket <j5j) . The noncanonical Euler-Lagrange equations (fTU| for the 
reduced gyrocenter coordinates (X,pii) are 



c dt X " dt u - v " sy ' 



X B * - b = Wfgy, (11) 

and 



b--S— = —SI = JL, (12) 
dt op\\ m 



where the gradient operator in Eq. (jlip is evaluated at constant pu and /i. 



2. Gyrocenter noncanonical Hamilton equations 

By inverting the Lagrange matrix w Q( 3 in Eq. (fT0|). we obtain the guiding-center Poisson matrix J a/3 = {Z a , Z^} sc 
from which Eq. ([5]) is constructed. Through this inversion, the gyrocenter Euler-Lagrange equations (|10[) are converted 
into the gyrocenter Hamilton equations 

dgyZ a _ -jafS d-Hgy _ 1 7 a tt \ _ ' ^ / R * tq/3 



where we have used the Liouville property of the guiding-center Poisson matrix J a @ and the guiding-center Jacobian 
Bp so that d a (Bt d gy Z a /dt) = 0. In addition, the gyrocenter noncanonical Hamilton equations for (X, py ) are 

rfgyP|| _ B* _ c ___ , _ / c _ c \ d gy X 
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•Vi/gy = -b-Vffgy + p„ (b-Vb) (15) 



The gyrocenter velocity (fT4|) contains the unperturbed guiding-center velocity p|| (B*/Sp -I- (cb/eBj|) X /iVB as 

well as the nonlinear perturbed gyrocenter E x B velocity (cb/BTA X V$ gy . The gyrocenter parallel force (fT5|). on 
the other hand, contains the unperturbed guiding-center parallel force — fi (B* / BT.) • V B as well as the nonlinear 
perturbed gyrocenter parallel force — e (B*/Bu) • V$ gy . 

Lastly, the gyrocenter magnetic moment y, is an invariant of the gyrocenter Hamiltonian dynamics, 



d gyA i _ n gggy 



0, (16) 



dt B d{ 

since the gyrocenter Hamiltonian ([2]) is independent of the gyrocenter gyroangle, while the gyrocenter gyrofrequency 

dgyC = fi af/gy = CO 9£gy 

dt B dfj, B dfj, { ' 

contains corrections to the unperturbed gyrofrequency Q that are due to the fluctuating electric field. 
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3. Gyrocenter canonical Hamilton equations 

In Sec. IIII1 we will need to express the gyrocenter Hamilton equations in canonical form which, with the help of 
the canonical gyrocenter momentum 

Pgy = + p,|b = - c A*, (18) 

satisfies the canonical Hamilton equation 



^gyPgy _ d-ffegy 
dt ~ <9X C [ 



(19) 



where the spatial gradient d/dX c is evaluated at constant p gy and (X, with H cgy (X c , /i, p gy ; t) = Hgy(X,p\i, (j,;t) and 
X c = X. Equation (fT9| is obtained as an Euler-Lagrange equation from the canonical one-form r cgy = p gy • dX c + 
Jd£ — H cgy dt, where the gyroaction coordinate J = (mc/e) /i is canonically conjugate to the gyroangle 

A direct application of the noncanonical Hamilton equations (fT4)) - (fl~5)) on the canonical gyrocenter momentum (fT8|) 
yields the Hamiltonian identity 

dgyPgy _ e d gy X ^ + ri gy p|| g _ £ / , ^ rf gy X + d gy X 




dt c dt dt c \ dt dt 

= -V*„ + fvA-.if S -^, (20, 

which relates the two different spatial gradients used in the Sec. IIIII [see Eq. (|94|)]. We henceforth use the notation 
(e/c) A* defined in Eq. (fT5)| whenever we use p gy at constant (e.g., dp gy = (e/c) dX • VA* + rfpn b). 

In Sec. IIV1 we will also need to write Eq. (|19[) in axisymmetric tokamak geometry, where the background magnetic 
field 

B = Vi^xVi/i + q(ip) V-0 x V?9 = V x A, (21) 

is represented in terms of the standard magnetic coordinates (ip, ip). Here, the safety factor q(ip) = B • V<p/(B • Vi?) 
is a function of the poloidal flux ip (note that B • Vip = 0), and the vector potential associated with Eq. (|2"TT) is 
A = — ip Vp + A$(ip) V#. From Eq. (fT8"|). the toroidal canonical gyrocenter momentum is therefore defined as 



9X _e 
gV ' Psy ~ ~c 



•Pgy = - - Ip + P\\ K> ( 22 ) 



where b v = b • dX/dp ~ 72. = |3X/oV| is approximately given by the major radius 72. of the axisymmetric tokamak 
plasma. The toroidal canonical gyrocenter momentum satisfies the canonical Hamilton equation 

^gy^ll h ( e dgyip ^syK\ _ d-ff gy <9<I> gy 



dt dt v \c dt P|1 dt ) ~ dp " dp ' 1 j 

where the last equality follows from the condition of axisymmetry on the background magnetic field. Equation (|23[) 
can also be obtained from the noncanonical Euler-Lagrange equation (jlll) after using the identity 

e dX d gy X _ e ( d sy tp d gy d 



X B* = -J I -2^- Vi? ^- J ■ B* 

c ay a£ c \ dt dt J 

_ £ '^gyV' _ d gy b v 
c dt " 

where J7" denotes the Jacobian for the transformation X — > (ip, $, <£>), which is defined from the identity 

J- 1 = (Vt/> X Vtf ) • V</? = B" 8 . (25) 

Equation follows from the identities 

dX dX dX dX 



^ OA OA ^ OA \ _ / ^ ^ \ 

dip ' OV di9 ) \ ) 
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with the magnetic-curvature terms 



™-»*.™.B-).(i-a$,a&) 



(26) 



which hold for the axisymmetric tokamak magnetic field (|2"Tj) . These identities are then used to construct d gy b v /dt = 
(db v /dip) d gy ifj/dt + {db v /dd) d gy ti/dt in Eq. (JUJ). 

Before moving on to the variational formulation of the gyrokinetic Vlasov-Poisson equations, we introduce the 
following gyrocenter equations, obtained from the noncanonical gyrocenter Hamilton equations (|14 j) -(|15 [) . that will 
appear in gyrokinetic Vlasov- moment equations presented in Sees. HITCl and [IV Al First, the gyrocenter equation for 
the parallel-momentum vector p|| b is 

d sy(P\\b) _ (B* \ r d gy b _ e d gy X d gy b 

.-V, gy + ^XB +P|| vi^, (27) 



c 



and, second, the gyrocenter equation for the toroidal component of the parallel-momentum vector pu b v is 

-dt {p ^ ] - e [c— t W) = I sy c^T xB J " v (28) 

where E gy = — V$ gy denotes the generalized gyrocenter electric field and we used the magnetic-coordinate identity 

<9X 

W = B x — , (29) 

satisfied by the axisymmetric tokamak magnetic field ([2~Tj) . 

Lastly, we introduce an important identity associated with the connection between the canonical gyrocenter mo- 
mentum Hamilton equation (|19j) and the canonical gyrocenter toroidal angular- momentum equation (|23|) . First, we 
take the dot product of Eq. (f2"U)l with dX/dtp and obtain 

^gyPgy . ^ X _ d gy p gyv e rf gy X /<9X . ^, 



dt c?f c dt \ dip 



c dip dt ' 
where the cylindrical dyadic tensor 

v(^) = 7^^ - (31) 

is defined in terms of the unit vectors 1Z = VIZ and (p = dX/dcp. Next, we insert Eq. (|23[) in Eq. (f30|) and obtain 
the gyrocenter canonical identity 



q _ d gy p gytp | dH gy _ e cigyX 



dt dip c dt 

By inserting Eq. pip into Eq. (|3"2"j) , with z = 1Z X <p, we obtain the identity 



dA* A 

dip \dip 



(32) 



which is trivially satisfied in the cylindrical representation A* = A^IZ + A* v (p + A*z, where the components are 
independent of the toroidal angle (p while the unit vectors?* = (1Z, (p,z) satisfy the identities de 1 /dip + e i Xz = 0. 
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B. Gyrokinetic Vlasov-Maxwell Variational Principle 

The derivation of exact conservation laws for the nonlinear gyrokinetic equations relies on the existence of variational 
principles from which they can be derived 0413 • We now briefly review the variational derivation of the gyrokinetic 
Vlasov-Poisson equations based on the principle of least action <5„4 gy = 0, where the action functional is [9j 

Azy = J ^ (VlExl 2 - |B| 2 ) - E / d s ZF gy (Z)H gy (Z-^) = J £ ey d*x, (34) 
where the gyrokinetic Lagrangian density C gy is the sum of the Maxwell Lagrangian density 

1 (e 2 |E X | 2 - |B| 2 ), (35) 
where the perturbed electric field is Ei = — V0i, and the gyrocenter Vlasov Lagrangian density 

£gyv = " £ / J sy( Z ) n ^ Z ' *P> (36) 

where the extended gyrocenter Hamiltonian 

n gy (Z; fa) = fl gy (X,p||, A1 ,t; 0i) - W 7 (37) 

is expressed in terms of the (regular) time-dependent gyrocenter Hamiltonian @ and the gyrocenter energy coordinate 
W. The extended phase-space integration in Eq. (|34|) is defined with d 8 Z = dtd 3 X d 4 p, where d 4 p = c _1 c£W d 3 p and 
d 3 p = TO-Bm dp|| diidC, (here, m_B| denotes the Jacobian). Furthermore, the extended gyrocenter Vlasov distribution 

T gy {Z) = c5(W-H ey ) F(X,p h »,t) (38) 
ensures that the gyrocenter Hamiltonian dynamics satisfies the physical constraint 7i gy = 0. 

1. Eulerian variations 

The gyrokinetic Vlasov-Poisson equations describing the self-consistent evolution of the gyrocenter Vlasov distri- 
bution F in the presence of electrostatic fluctuations </>i are obtained from the Eulerian variational principle 

SC gy d 4 x = 0, (39) 

where the Eulerian variation of the gyrokinetic Lagrangian density is expressed as 



5C gy = £ (flVEx) -Y.J 



d*p, (40) 



with d 3 p = 2irmB^ dp\\dfi (which henceforth assumes gyroangle-averaging) . In addition, B is not a variational 
field (since it is time independent) while the Eulerian variations <5Ei = — VStfii preserves the electrostatic constraint 
V X 5~Ei = 0. The Eulerian variation for the extended gyrocenter distribution ([38]) is 5T gy = {S gy , J- gy } gc , which 
preserves the Vlasov constraint J 8F gy d 8 Z = under a virtual canonical transformation Z a — > Z a + SZ a in extended 
phase space, where SZ a = {Z a , S gy } gc is generated by the extended scalar field S gy . 
Equations ([2]) and (|SJ) yield the functional derivative 



<tyi(x) 



(T-, 1 5 3 gc ) = e (,5 g 3 c - e {S lt S 3 gc } gc + •■•), (41) 



where the guiding-center delta function <5 gc = 5 3 (X + p gc — x) indicates that the gyrocenter contribution at a fixed 
field point x only comes from gyrocenters located on the ring X = x — p gc . It is convenient to formally express the 
combined guiding-center and gyrocenter push-forward operations 

T" 1 ^ ee T^fX - x) ee 5 3 (X + p e - x) (42) 
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in terms of the reduced displacement 



Po 



Pi 



(43) 



where the unperturbed gyrocenter displacement p = p gc0 denotes the (lowest-order) guiding-center gyroradius (which 
is gyroangle-dependent) and the first-order reduced displacement 

Pi = esP gc i + ePg y i = (Pi) + Pi (44) 
includes the first-order guiding-center displacement p gcl (22l . [23l | and the first-order gyrocenter displacement 



Pgyl = \ X + p g( ., Si \ , 
I. J g C 



(45) 



which is expressed exclusively in terms of 5*1 (and therefore 0i gc ) according to Eq. ([8]). Each first-order contribution 
in Eq. (|44j) has a gyroangle-independent part, which contributes an electric-dipole term in the gyrokinetic polarization 
[see Eqs. (|54[) -(j57 p ], and a gyroangle-dependent part, which contributes an electric-quadrupole term in the gyrokinetic 
polarization [see Eq. (|55)) ]. (We note that effects due to the gyrokinetic magnetization [10] do not enter into our 
present discussion but will be considered in future work.) 

By rearranging terms in the Eulerian variation (|4T)|) so that the variation generators (iS gy , 5(j>i) appear isolated, we 
obtain 



e 

47T 



E 



V-Ei - e J Fgy (Tgy 1 ^) d 4 pd 3 X 
Sgy { 



' -^gy ' ^gy r ^ P 



<9A „ „ 



where the Noether fields A and T are 

A -12 



^gy -^gy d P> 



e 

47T 



k Ei 



^£ry -^gy 



at 



(46) 

(47) 
(48) 



We note that the Noether space-time divergence terms dA/dt + V-T do not contribute to the Eulerian variational 
principle but instead play a crucial role in the derivation of exact conservation laws by the Noether method (see 
SecHIIl. 



2. Gyrokinetic Vlasov-Poisson equations 

When the variation (1461) is inserted into the Eulerian variational principle (1391) for arbitrary variation generators 
(<S gy , 54>i), we obtain the gyrokinetic Vlasov equation in extended phase space 

{Jgy, Hgyjge = 0, (49) 

and the gyrokinetic Poisson equation 

eV-Ei = 4tt^ e J F gy (T" 1 8%) d 4 pd 3 X. (50) 

Furthermore, when the energy integration ( J dW) is performed on the extended gyrokinetic Vlasov equation (|49[) . we 
obtain the gyrokinetic Vlasov equation 

dF r i dF d EV X eLyPii dF , x 

9t I ' sy / g c <9t eft di ^ ; 

where dF/ dfi is absent because of Eq. (TT6|) and the gyrocenter gyrofrequency (pTT|) does not appear in Eq. (fBTj) since 
dF/dQ = (by construction). Hence, the guiding-center and gyrocenter transformations yield a reduced dynamical 
Hamiltonian description of gyrocenter motion in terms of four-dimensional coordinates X and pii , with fi appearing 
as a parameter. 
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The gyrokinetic Poisson equation (|50[) . using Ei = — V0i, becomes 
where the gyrocenter charge density is defined as 



6gy 



E 



F<Tp, 



and the gyrokinetic polarization 



(52) 



(53) 



(54) 



includes the electric-dipole contribution associated with the first-order reduced displacement (|43[) . where (p £ 
(pj) + • • • (since (p ) =0), as well as a quadrupole contribution represented by the second-rank tensor 



Q-E 



(55) 



which includes standard finite-Larmor-radius corrections due to (PoPo) as we ^ as the first-order corrections (p Pi) 
and (p 1 p ) associated with the gyroangle-dependent part of the first-order displacement (|4"4")l . The first term in the 
gyrokinetic polarization (|54j) includes the first-order guiding-center contribution (2(| 



- n x dt 

and the first-order gyrocenter contribution 

(Pgyi) = 



^ V±B H !!■ b- Vb 

m 



e d I I 
B djl \ Po ^ 



By using Eqs. ()56P -([57 |) . the gyrokinetic polarization (|54|) can thus be expressed in the form 



do 



F 



dt 



d 3 p 



V-Q, 



where the truncated gyrocenter velocity 



dt 



eB 



e B \ [i VB + — b • Vb 



eeV( 



(56) 



(57) 



(58) 



(59) 



denotes the perpendicular components of the gyrocenter velocity (fl"4"|) with the effective gyrocenter potential (U]) 
replaced with its first-order contribution e (4>i gc )- 



III. GYROKINETIC MOMENTUM CONSERVATION LAWS 

In this Section, we use the Noether method to derive an exact momentum conservation law for the gyrokinetic 

Vlasov-Poisson equations ([5"l"j) - ([52"T) . After substituting the gyrokinetic Vlasov-Poisson equations (|49|) -([50 |) into the 
Eulerian variational equation f|46[) , we obtain the gyrokinetic Noether equation 

6£ gy = ^- + V-T, (60) 

where the variations (iS gy , Sipi ) in Eqs. (|47l) - (|48[) are no longer considered arbitrary but are instead generated by 
infinitesimal space-time translations or rotations. The energy conservation law for the gyrokinetic Vlasov-Maxwell 
equations was obtained from the gyrokinetic Noether equation (jBT))) and was discussed in Refs. (3-fTolj by considering 
infinitesimal time translations t — > t + St. 



9 



Noether Method 



The momentum conservation law for the gyrokinetic Vlasov-Poisson equations (|51l) - ([52|) is obtained from the gy- 
rokinetic Noether equation (|60p by considering infinitesimal space translations x — >• x + <5x [H, |7l|, [27j , for which 



Pgv • Sx. 



gy — fgy 

5<f) 1 = -<5x-V(/>i = 5x-Ei 
SC ey = - V- (<5x£ gy ) + Sx-VCgy J 



(61) 



where the gyrocenter generating scalar field <S gy generates a virtual spatial translation (Sx and the gyrocenter canonical 
momentum p gy is defined in Eq. (fT8|) . 

The first term in the expression for 5C gy in Eq. (|61[) takes into account the geometric nature of £ gy as a space-time 
density. The notation V'£ gy in the second term, on the other hand, denotes the explicit spatial gradient of the 
Lagrangian density £ gy with the gyrokinetic fields (F, Ei) held constant (i.e., V = in a uniform magnetic field). 
Hence, the explicit gradient V'£ gy becomes 



B 

4tt 



E / F (v'h s 



(62) 



where the first term represents the contribution from the Maxwell part Q35[) of the Lagrangian density while the 
gyrocenter Vlasov term involves the explicit canonical gradient expressed in terms of the Hamiltonian identity (|20[) . 

By inserting the variations (|61[) into the gyrokinetic Noether equation (|60[) . we obtain the primitive gyrokinetic 
Noether momentum equation 



fa-V'£ gy = | 



E 



F ( Pgy • H d*p 



E / (Pgy • 5x ) <* 3 p - ^ (<**• Ei) Ei + 5x 



(63) 



where only the Maxwell part (|3"5j) of the gyrokinetic Lagrangian defined in Eq. (I3~4l survives the physical constraint 
■H gy = 0. Next, if we replace Cm and V'£ gy in Eq. (1531 with Eqs. (1331) and (|6"2")l , respectively, we obtain 



- <5> 



E 



c dt 1 



d_ 
dt 



where the right side involves the gyrocenter canonical momentum p gy and the perturbed electric field Ei. 



(64) 



B. Gyrokinetic Linear Momentum 

The gyrokinetic Noether (linear) momentum equation is obtained from Eq. (|64[) by considering a virtual translation 
generated by an arbitrary constant displacement 5x. The gyrokinetic Noether momentum equation is, therefore, 
expressed as 



<9P 



where the gyrokinetic momentum density is 



^ + v.n = -£/F(v'H, 



, .£vA-.*£)*ft 



E / F Pgy d3 P> 



(65) 



(66) 
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and the canonical gyrokinetic momentum-stress tensor is 

n = i ( |Ei|a \ EiEi ) + £ J F %r Pgy d3p = He + Hgy (67) 

Here, the canonical stress tensor contains the symmetric Maxwell stress tensor 

n E = ±- (|Ei| 2 I - Ei Ei) , (68) 



and the asymmetric gyrocenter canonical momentum-stress tensor 

n - - £ / F p - dV (69) 

We note that, while the asymmetry of the gyrocenter canonical momentum-stress tensor (|69p is not physically relevant 
here since it appears in the gyrokinetic Vlasov- moment equation (j78|) , which has no symmetry requirements on n gy , it 
plays a crucial role in deriving the exact gyrokinetic toroidal angular-momentum conservation law [see Eqs. (|92[) - (|94[) ]. 

We now show that Eq. (|55|) is consistent with Eq. (ITSl by constructing an explicit proof of the gyrokinetic Noether 
momentum equation (|65[) . We begin with the partial time derivative of the gyrokinetic momentum density (|66[) 



dt ^ J \ dt dt 
Upon substituting the phase-space divergence form of the gyrokinetic Vlasov equation (fBTj) . Eq. (fTOf becomes 

Next, we compute the divergence of the canonical momentum-stress tensor (|67p to obtain 

V-n = -e ^ e J F (TgyEigc) d 3 p + V-II gy , (72) 



where we used the gyrokinetic Poisson equation (|50|) with Ei = — V0i and Ei gc = — V<^i gc , and we used the 
Maxwell-tensor divergence 

V • n E = - e e / F ( T gy lE igc) rf 3 P- (73) 
By inserting Eqs. (|7i p ~(|72 p into Eq. (1651) . we obtain the gyrokinetic momentum equation 

d 3 p. (74) 



'^ P - ee <T^E lgc ) + V'tf gy - - VA* X 



dt \-gy— sw C dt 

We immediately conclude that the canonical gyrocenter momentum p gy must satisfy the gyrocenter equation of motion 

d g.vPgy _ r „ /t-Itj. \ (u'u e \ * d sy X 



A -e(T-Ei gc >-^ gy --VA*.^j. (75) 
We note, however, that the expression for (TgyEigc) can be rewritten as 

e (Tg-^Exge) = -eV(0 lgc ) + e 2 V0 lgc }) + ••• 

V$ gy - V'$ gy ) , (76) 



where the explicit gradient term V'<I>g y takes into account the nonuniformity of the background magnetic field (e.g., 
which appears in the guiding-center Poisson bracket). Lastly, using the relation eV'f gy — V'H gy = — Vi? gc , we 
immediately recover the canonical gyrocenter momentum equation (|19p from Eq. (|T5|) : 



dt -"I V sy c 

^ = - V/, u , 

dt &y c dt 



e „ . j, <L V X\ „„ e . „ cL v X 



= - e eV$ gy - (VHg C - :VA».^) = - VH gy + -VA*-^p. (77) 
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We therefore see that the gyrokinetic Noether momentum equation (|65|) is automatically guaranteed by the canonical 
gyrocenter momentum equation (fT9"|. Lastly, we note that Eq. ([71)) can be written as 



^gyPgy 
dt 



e _ . ± cL y X 

- VA* • 

c dt 



d 3 P 



§ + v .n gy + i:/F(vH g) .-£vA..^)^ 



which can be recovered as a gyrocenter Vlasov-moment equation 



9 

°= at 



xFd P 



V- 



at 



obtained from the gyrokinetic Vlasov identity 







d_ 

at 



(B; X F) + —^ X F{ Z ",H gy }J-BlF 



dgyX 
dt ' 



(78) 



(79) 



(80) 



which holds for any function \ on gyrocenter phase space (e.g., \ = Pgy)- We note that Eq. (|80p follows from writing 
the gyrokinetic Vlasov equation (|51|) in phase-space-divergence form using Eq. (fT3|) . 



C. Gyrokinetic Parallel Momentum 

As a simple application of the gyrokinetic Noether momentum equation (|65j) . we now derive the gyrokinetic parallel 
momentum equation by first decomposing the gyrokinetic momentum density (|66[) and the gyrocenter momentum- 
stress tensor (|6T?1) as 

(P, n gy ) = ( p,| + - c ggy a, n gy ,| + - c J gy a) , (81) 

where g gy is the gyrocenter charge density (|53l) and the gyrocenter current density is 

J g y - E e / F ^ d ^ ( 82 ) 
while the gyrocenter parallel momentum vector Pm and the gyrocenter parallel-momentum stress tensor II|| are 

(P||,n gy|| ) = £ | F (S, ^b) P|| d 3 P . (83) 

Secondly, we introduce the decomposition 

e _ . * cL y X e d CT X cL v X 
- - VA* • = VA • -Sf pi, Vb • - §y — 

c dt c dt dt 

and the Maxwell-tensor divergence 

V-II E = -e £ e J F (T^Exge) d 3 p = £ e | F (V$ gy - V'$ gy ) d 3 p, 
which combines Eqs. (|73|) and ()76|) . Thirdly, we use the gyrocenter charge conservation law 



= "# + v * Jgy = dt \ Sphys + v ' v ) + v ' \ Jphys ~ ~dt ~ cVxM 
= ^dT + v ' Jphys (84) 

where the physical charge density g p h ys = (?gy — V • "P is expressed in terms of the gyrocenter charge density (|53[) and the 
gyrokinetic polarization charge density — V • "P, while the physical current density J p h ys = Jgy + dP /dt + c V X M is 
expressed in terms of the gyrocenter current density J gy , the polarization current density dP/dt, and the divergenceless 
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magnetization current density c V X M (note that the gyrokinetic magnetization M is not needed in the present 
gyrokinetic Vlasov-Poisson theory). Hence, we obtain 



do. 



gy 



A + V.(j gy A) = J gy .VA = £e / F^-VA. 



dt 

By inserting these expressions, and using Eq. (|27[) . we obtain the vector gyrokinetic parallel- momentum equation 

— + V.n gy|| = ^ / F^pdp, (85) 



which is clearly expressed as a gyrokinetic Vlasov-moment equation (|T9[) with x = b. 

Lastly, we take the dot product of Eq. (|85l) with the parallel unit vector b and we obtain the gyrokinetic parallel 
momentum equation 

^ + v.rii = - E / - Pn (s.vb) .3s*) d^p^Y, J F ^ d ^ ( 86 ) 

where (P||,R||) = (P|| • b, ITg y || -b) and we used (V-II gy ||) -b = V-R|| with Vb : IlJ y || = (since Vb-b = 0). We 
immediately note that the gyrokinetic parallel momentum equation (|86[) has the form of a gyrokinetic Vlasov-moment 
equation (|79|) with \ — P\\ ■ 



D. Gyrokinetic Toroidal Angular Momentum 

We are now ready to use the Noether method to derive the gyrokinetic toroidal angular-momentum conservation 
law associated with the toroidal rotational symmetry of the background magnetic field. First, we insert a virtual 
infinitesimal toroidal rotation (about the z-axis) generated by the displacement 

<5x = — Sip = 6(pz X X (87) 
dip 

into the primitive gyrokinetic Noether momentum equation (|63[) . The gyrokinetic angular- momentum conservation 
law is, therefore, expressed as 

BP 

+ v-n v = 0, (88) 
where we used the axisymmetry condition d'£ gy /d<p = 0. Here, the gyrokinetic angular- momentum density is 

p * = E / d3 p> ( 89 ) 

where the gyrocenter toroidal canonical momentum p gyip in defined in Eq. (|22f) , and the canonical gyrokinetic angular- 
momentum flux 



<9X 

n v = n E — + n v (90) 



is decomposed in terms of the toroidal component of the Maxwell stress tensor Ip5|) and the toroidal component of 
the gyrocenter Reynolds stress tensor 

R ^ E/ F^P S y,d*p, (91) 

where the gyrocenter velocity is defined in Eq. (|T4|) . We note that the connection between the gyrokinetic angular- 
momentum conservation law (|88p and the gyrokinetic Noether momentum equation (|65p involves the dot product of 
the gyrokinetic Noether momentum equation (|65p with d~K/d(p: 

f-f + <v.n,.f ^ + v.n,-v(f):n T 
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where II T denotes the transpose of the canonical gyrokinetic momentum-stress tensor (j67]) . Noting that the dyadic 
tensor (1311) is antisymmetric, only the asymmetric part of the canonical gyrokinetic momentum-stress tensor (|67[) 
contributes to the term 



ax 



n 



v 



ax 

dp 



n 



E 



F 



e d gy X 



dt 



dip 



•A* 



d 3 P , 



(93) 



since the Maxwell stress tensor IIe is symmetric and V(<9X/<9<^) : II J = 0. By substituting Eq. (|93[) into Eq. (|92[) . 
we recover the gyrokinetic toroidal angular-momentum conservation law 







dP v 
dt 



vn 



E 



F 



e dgyX 



<9A* 

dp> 



dX 

dp 



where the gyrocenter canonical identity (f3"2"j) reappears. 

Lastly, we note that the toroidal component of the Maxwell-tensor divergence (|73|) yields 



(v-n £ 



gx _ Y> 



„ 9$gy ,3 

F _gy d 3 pj 
d<p 



where 



d 3 p, 



(94) 



(95) 



dp 



<9</J 



5i, 



3lgc 



dp 



T 



c)0 



lgc 



dp 



When we combine these contributions into the gyrokinetic angular-momentum conservation law 



dP v 
dt 



E 



F 



/dH g 



V dp 



d 3 P , 



(96) 

we finally obtain 
(97) 



where dH gy /d = e d§ gy /dp. Note that Eq. (j^T)) can be obtained from the gyrokinetic Vlasov-moment equation (|T9")) 
with x = P eyi p-- 







E 



F 



gy-Pgyy 



dp 



d 3 p 



dP v 
dt 



V-R, 



E 



F 



/ml 



V dp 



d 3 p, 



as was recently derived by Scott and Smirnov [l7|. Equation (l97|) is the general form of the gyrokinetic angular- 
momentum conservation law and applies to the electromagnetic case as well as the electrostatic case considered here. 
In the electromagnetic case (to be considered in future work), Eq. (|97|) applies to the Hamiltonian formulation of 
electromagnetic gyrokinetic theory [H, Q , where all perturbation fields appear in the gyrocenter Hamiltonian and the 
gyrocenter canonical momentum is still given by Eq. (j 18[) . In the next Section, we will investigate this conservation 
law in axisymmetric tokamak geometry (|21j) for the case of the gyrokinetic Vlasov-Poisson equations. 



IV. GYROKINETIC MOMENTUM CONSERVATION IN AXISYMMETRIC TOKAMAK GEOMETRY 



In this Section, we derive the gyrokinetic toroidal angular-momentum conservation law (|97p in axisymmetric toka- 
mak geometry (|2ip . For this purpose, we introduce the magnetic-surface average 

[■••I = ^ <f (■■■)Jd$dp, (99) 

where V(ip) = § J dd dpi is the surface integral of the magnetic-coordinate Jacobian. The flux-surface average (|99|) 
satisfies the property 

(100) 

for any vector field C. In a time- independent axisymmetric tokamak geometry, we note that d/dt also commutes with 
magnetic-surface averaging. 
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A. Gyrokinetic parallel-toroidal momentum transport equation 



We now present a derivation of the magnetic-surface-averaged gyrokinetic parallel-toroidal transport equation (|97l 
in axisymmetric tokamak geometry, carried out in four steps. 



1. Parallel-toroidal canonical momentum 



As a first step, we begin with the gyrocenter toroidal canonical momentum (|22j) . where we quickly note here that 
only the parallel component of the gyrocenter toroidal angular velocity dip/dt appears (i.e., the magnetic-drift velocity 
is absent from Pgy V ). We therefore write the gyrokinetic angular- momentum density (|89|) as 



R, 



d 3 p 



£gy + P \\v-> 



(101) 



where g gy denotes the gyrocenter charge density (|53l) and the parallel-toroidal gyrocenter momentum density is 



Pii 



J Fp {] d 3 p\ b, 



The partial time derivative of Eq. (|101[) yields 



tp dg 



gy 



dP, 



Ik 



dt 



so that Eq. (j9"7| becomes 



dR 



dt 



c <9i 
c 0t 



E 



<9Hgy 3 

F - sy d^p, 
dip 



(102) 



(103) 



(104) 



which represents the evolution equation for the gyrocenter parallel-toroidal momentum density P\\ v , where H v is 
defined in Eq. (|9T]). 



2. Magnetic- surf ace average 



As a second step, we perform the flux-surface average of Eq. (|104|) . which yields 



dt 



1 d_ 

V d^ 



(v Kl) 



d[ggyl 

c at 



E 



F^hL^p 
dip 



(105) 



where the surface-averaged radial flux density B?{t = • Vf/ 1 of gyrocenter toroidal canonical momentum 



E 



F 



eft 



p gyip d p 



is decomposed in terms of the radial component J^, = J gy • Vt/> of the gyrocenter current density 
flux of gyrocenter parallel-toroidal momentum 



E 



dt 



By rearranging terms in Eq. (|105[) . we therefore obtain 

d[P\\vl 1 ± 
dip 



dt 



V dib\ 



Rf 

\\v 



= E 
= E 



PH (Pp 6 



. 1 dgyV 



9y> 



d 3 p 



F 



dt 



(Pll M 



d 3 p 



where we used the surface-averaged gyrocenter charge conservation law 

dlggy] - _ t v .j 1 = _IA(v 



(106) 
and the radial 

(107) 



(108) 



(109) 



By using Eq. (|2"5|) . we see that Eq. (|108l) is expressed directly as a gyrocenter Vlasov-moment equation (fT9")l with 

X=P\\ V 
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3. Gyrocenter quasmeutrality condition 



As a third step, we make use of gyrocenter quasineutrality relation 

£> ev = V-V 



(110) 



between the gyrocenter charge density (|53|) and the gyrokinetic polarization (|54[) . The partial time derivative of the 
surface-averaged gyrocenter quasineutrality condition (|110j) yields 



Of 



dV 
~dt 



dt 



an) 



where = "P • Vip denotes the radial part of the gyrokinetic polarization and dV^ /dt defines the radial component 
of the gyrocenter polarization-drift current. 

By using the surfaced-averaged gyrocenter charge conservation law (|109l) , we readily obtain the identity 



dt 



0. 



(112) 



which implies that the surfaced-averaged radial component of the physical current JJ p hy S ■ V^] = vanishes, as 
is required to preserve the ambipolarity condition [l9l |28| . Here, we used the fact that the surface-averaged ra- 
dial component of the magnetization current density c V X M vanishes since V-0 • V X M = V • (M X V^) and 
JV • (M X W)] = by Eq. (fIDD)> . 



4- Parallel-toroidal momentum transport equation 



As a fourth and final step, we substitute Eqs. (|109[) and (1112[) in Eq. (|105l) . and we obtain the gyrokinetic transport 
equation for the surface-averaged gyrocenter parallel-toroidal momentum density 



d_ 

dt 



ld_ 

V£hp 



E 



dip 



(113) 



which corresponds exactly to Eq. (80) of the recent work by Scott and Smirnov [17| . Here, we see that the transport 
of the surface-averaged gyrocenter parallel-toroidal momentum density [-Piu] is combined with the surface-averaged 
radial gyrokinetic polarization \P^\. In fact, using the magnetic-coordinate identity (|29|) and the gyrokinetic polar- 
ization (158p . we find that the sum of the gyrocenter parallel-toroidal momentum density and the radial component of 
the gyrokinetic polarization 



F 



'ax 

dp dt 



d 3 p 



V^- (V-Q) 



(114) 



is expressed in terms of the gyrocenter moment of the toroidal component of the gyrocenter momentum (with the 
effective gyrocenter potential $ gy replaced with its first-order contribution e (4>i gc }) as well as quadrupolc contributions 
to the gyrokinetic polarization. The combination of the gyrocenter canonical toroidal momentum equation (1231) and 
the gyrocenter charge conservation law (|111[) . therefore, allows us to recover the gyrocenter Vlasov moment equation 
for the (truncated) toroidal momentum 



^gy X 
dt 



dX 

dp 



— = mil* 



d^p 
dt 



(115) 



We note here that previous works (e.g., Ref. [17|) have neglected the guiding-center polarization contribution in 
Eq. (|58p and, hence, the first term on the right side of Eq. (|114l) only combines the toroidal components of the 
gyrocenter parallel vector momentum and the perturbed E x B velocity and ignores the toroidal component of the 
magnetic drift velocity. 
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B. Gyrokinetic parallel-toroidal momentum conservation law 



In their recent paper [17j . Scott and Smirnov considered a simplified gyrocenter Hamiltonian obtained in the long- 
wavelength (zero-Larmor-radius) limit, where the gyrocenter potential Q became $ gy ~ e<f>\ + \ e 2 p 1± • Vj_</>i and 
the first-order gyrocenter displacement (1571) was {p\, rY ) — P\±_ = — (c/BCl) Vj_0i. In this limit, Scott and Smirnov 
flTl found 



dtp 



dip 



e 2 Pu 



■ V 



d<p 



and then showed that the last term on the right of Eq. (|113[) can be expressed as a spatial divergence [as can already 
be seen from Eq. (|95p]. 

In the present paper, we perform a Taylor expansion of Eq. (|96[) in powers of the gyrocenter displacement p el 
defined in Eq. (pt5)) . which yields 



dH t 



dip 



sy = ee 



T _i (dfa 
6 {dip 



[dip 



(Pe)-V 



dip 



(116) 



where we have combined the guiding-center and gyrocenter phase-space transformations T e 1 = TgyTg,, 1 [see Eq. (|42p] 
From this expression, we now obtain 



E 



F 



9H^ 

dip 



dp = eg, 



(117) 



Next, we substitute the definition of the gyrokinetic polarization (|51|) into Eq. (|117[) and, after rearranging terms, we 
obtain 



E 



<9-ffgy 3 

F „ sy cTp 
dip 



eV 



<9</? 



90i ^ 
dip 

eQ- V 



(118) 



where we used the gyrocenter quasineutrality condition (jl 10[) to show that the last term on the right of Eq. (|113[) can 
indeed be expressed as an exact spatial divergence involving moments of the gyrocenter displacement p e . We note 
that the right side of Eq. (|118[) , which is at least of first order in e, retains all finite-Larmor-radius effects as well as 
nonlinear corrections, such as the gyroangle- independent part (|5T[) of the first-order gyrocenter displacement Pi gy , 
while Scott and Smirnov [l7| considered only the long- wavelength limit of Eq. (|118[) , with the gyrocenter contribution 
to the gyrokinetic polarization T ~ — e (mnc 2 / B 2 )X7 ±<pi defined in the zero-Larmor-radius limit and all higher-order 
effects were omitted [see Eqs. (59) and (91) of Ref. [17J]. 

When we insert Eq. (|118[) into Eq. (|113l) . we obtain the gyrokinetic parallel-toroidal momentum conservation law 
for the nonlinear gyrokinetic Vlasov-Poisson equations: 



- IP*} 

c 



ld_ 



dip dip 



= 0, 



where the radial gyrocenter quadrupole vector 



J2lJ F{ PePe )d^ = Qt 



+ 



(119) 



(120) 



includes guiding-center FLR corrections Qq (involving p ) as well as higher-order gyrocenter corrections Q{ + • • • 
that involve the gyroangle-dependent part of the first-order guiding-center displacement p gcl and the first-order 
gyrocenter displacement p gy i (e.g., through (p Pi))- The gyrokinetic parallel-toroidal momentum conservation law 
(|119p therefore implies that the transport of gyrocenter parallel-toroidal momentum [-Pii^] is intimately connected 
to the gyrocenter radial polarization fPj\ . This connection has recently been investigated through the gyrokinetic 
Vlasov- moment approach [l2l [Til [TBI . Il7| . The term d(j>i / dip} in Eq. (|119l) represents the turbulent transport of 
toroidal momentum [l2l.ll3|. while the term [Q - ^ • Vd(j>\/dip\ involves the radial gradient in the toroidal electric field. 
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TABLE I: Gyrokinetic Vlasov-moment Equation (I79[) for various functions \- 



Function \ 


Pgy 


Psy<p 


PI 


b 




pw 


Gyrokinetic Vlasov-moment Equation 


Eq. (EH) 


Eq. J98J 


Eq. 


(ESI 


Eq. (nnu 


Eq. ((861) 



Lastly, we note that Eq. (|1 contains the following nonlinear ordering 

— 1 T c 
e w fi • (mc s U)P\\ ~ - • e(k±p s )ep s ■ (k ± p s )—, (121) 
a eq 

where, on the left side of Eq. (|121[) . O -1 d/dt = e u defines a turbulent time-scale ordering parameter and [P^] ~ 
(mc s 7V) Pj| defines the normalized parallel momentum, while, on the right side, |V^| d/dip ~ 1/a = e/p s introduces 
the ordering parameter e = p s /a, P^/|VV> ~ e (k± p s ) e p s represents the ordering of the gyrocenter contribution to the 
radial gyrokinetic polarization (|57|) . and ed<p\/d<p ~ (k±p s /q)T e /e represents the ordering of the toroidal derivative 
of the perturbed electrostatic potential (here, we used e0i ~ eT c /e, defined in terms of the electron temperature 
T e , and d/dip ~ k±a/q = k±p s /qe, where we assume fcj^ a 3> fcii 7?.). Hence, we obtain from Eq. (|12ip the following 
dimensionless ordering 

e w P|, ~ e 3 (fc ±(0s ) 2 a /(^). (122) 
If the turbulent time-scale ordering e u satisfies a gyro-Bohm (gB) ordering (2t| 



then the normalized parallel momentum P|| satisfies the explicit scaling 

P,l ~ (k ±Ps ) 2 a/(q1Z), (124) 

which involves the fluctuation spectrum in k± and physical parameters of the axisymmetric tokamak plasma. For 
typical modern tokamak plasmas, with c s ~ 1000 km/s (at T c ~ 10 keV) and a/(qlZ) ~ 1/10, we can insert the 
scaling (|124[) into the experimental upper-limit on intrinsic rotation speeds U\\ ~ Pii c s < 40 km/s and obtain 
the condition (/c_l/9 s ) 2 < 0-4 on the fluctuation spectrum, which is consistent with the standard gyrokinetic ordering 

. IK < 1 1- 



V. SUMMARY AND DISCUSSION 



In the present work, we presented the derivation of the exact gyrokinetic Noether momentum equation (|65[) and 
the exact gyrokinetic toroidal angular-momentum conservation law (|88[) associated with the nonlinear gyrokinetic 
Vlasov-Poisson equations. While previous gyrokinetic momentum transport equations were derived from moments 
of the gyrokinetic Vlasov equation, the gyrokinetic momentum equations derived here are exact statements which 
depend on the nonlinear gyrokinetic physics included in the gyrokinetic action functional (|34|) . Various proofs were also 
presented throughout the paper to demonstrate the interconnections between the gyrokinetic momentum conservation 
laws (p?2"j) and (|94p , and several intermediate steps were shown to be expressed as gyrokinetic Vlasov-moment equations 
dZnj) (see Table©. 

The gyrokinetic parallel-toroidal momentum conservation law (|119p , derived in axisymmetric tokamak geometry and 
averaged over magnetic-flux surfaces, highlights the role of polarization effects that appear explicitly in the surface- 
averaged gyrokinetic radial polarization JP^] (which includes guiding-center and gyrocenter contributions) and the 
toroidally-dependent parts that appear in [P^ d(f>i / dip + • Wd(j>i /dip +■•■]. The major difference between previous 
works (e.g., Ref. [ij]]) and our work is that the gyrokinetic polarization (|55|) considered here includes contributions 
from the guiding-center and gyrocenter phase-space transformations, while previous works have neglected the guiding- 
center polarization contribution (which is of the same order as the gyrocenter contribution). 

In future work, the exact gyrokinetic momentum and angular momentum conservation laws for the nonlinear 
gyrokinetic Vlasov-Maxwell equations, where fully electromagnetic fluctuations as well as gyrokinetic polarization 
and magnetization are taken into account, will be derived by the Noether method. Because the gyrokinetic Noether 
momentum equation (|65p was derived in general magnetic geometry, applications of the present work to stellarator 
plasmas [19h21| and other nonaxisymmetric magnetic geometries could also be considered. 
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Appendix A: Noether Method 



In this brief Appendix, we use the Noether method to derive the energy-momentum and angular-momentum con- 
servation laws [3lTl33l ]. In order to simplify the presentation, we assume the Lagrangian density C(ip a , <9 M -0 a ; x, t) 
depends on a iV-component field ijj a = (ip 1 , ip N ) and its space-time derivatives d^ip 11 = (c~ 1 dt4> a , Vip a ) as well as 
possible space-time dependence due to external fields. Furthermore, we assume that each component ip a , which is to 
be varied independently of each other, satisfies the Euler-Lagrange equation 



d ( dC 



dC 



dx^ \d{d fJi il) a 

In this case, the Noether method is based on the Noether equation 
dC dSip a dC 



(Al) 



SC 



5ip a 
d 



dip a dx^ d{d^ a ) 



S^ a 



' dC 


-A- 


( dc V 


d 


' dC 




dx^ 









Stp a 



dC 



S^ a 



(A2) 



where the variations 8ip a and SC will now be expressed in terms of infinitesimal space-time translations (energy- 
momentum) or infinitesimal rotations (angular- momentum). 



1. Energy-momentum Conservation Law 



In order to derive the energy-momentum conservation law through the Noether method, we consider arbitrary 
infinitesimal space-time translations generated by Sx v : 

x v -> x v + Sx u . (A3) 

In this case, the variations 5tp a and SC are expressed as 

5xj; a = -Sx v d u ijj a , (A4) 
SC = — d u (Sx v C] + Sx" d' v C, (A5) 



where the partial derivative d' v C is calculated with the fields tp a held constant. If we substitute the variations ()A4|) - (|A5j) 
into the Noether equation (|A2|) , we obtain the energy- momentum equation [33| 



<9 M T^ = d'"C, (A6) 
where the canonical energy-momentum tensor T^ v is defined as 

dC 

T v = g v c - — — - d v r- (A7) 
o{d l _ L %p a ) 

The Noether Theorem now states that if the Lagrangian density C is invariant under space-time translations along 
the x x coordinate (i.e., d' x C = 0), then the components T^ A satisfy the energy-momentum conservation law 

= d^ x = c" 1 d t T 0X + djT jX . (A8) 
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For time-translation invariance (A = 0), we define the energy density £ = T 00 and the energy-density flux = cT J °, 
and Eq. (|A8p yields the energy conservation law 

^+V-S = 0. (A9) 
at 

For space-translation invariance (A = j), we define the momentum density P J = T 0j /c and the stress tensor T y = IP- 5 , 
and Eq. (|A8I) yields the momentum conservation law 

<9P 

— + v-n = o. (Aio) 

2. Angular-momentum Conservation Law 

In order to derive the angular-momentum conservation law through the Noether method, we consider arbitrary 
infinitesimal rotations generated by the antisymmetric tensor 5A ua = — 5A av : 

x v -¥ x v + SA ua x a . (All) 

In this case, the variations 5ip a and SC are expressed as 

6^ a = - i 5A ua (x a d^ a - x v d a i> a ) + A^ a , (A12) 

5C = - d u (5A ua x a C^j + i SA"" (x a dlC - x v d' a Cj . (A13) 

In Eq. (|A12|) , the Lagrangian variation Aip a associated with the infinitesimal rotation (jAllj) is defined as 

Ar = ~ SA^ lj« u , (A14) 

where the antisymmetric tensor (whose role will be clarified below) depends on the field components (ip 1 , ip N )- 
If we substitute the variations (|A12I) - (|A13|) into the Noether equation (|A2[) . we obtain the angular-momentum 
equation 

q^jhI""] = x °d w L - x"d' a £, (A15) 
where the total angular-momentum tensor J^ 1 " 7 ^ = — J^\- au \ is defined as 

jli[ua] = J^lAf"?] _|_ g^[v<r}^ (A16) 

Here, the orbital angular- momentum tensor 

L* ff l ee x a - x v , (A17) 
where the canonical energy- momentum tensor T^ v is defined in Eq. (| A7|) . satisfies the equation 

The spin angular-momentum tensor 



S^ [Va] = 7777< ~LO a ^ (A19) 



a^S^l ee T va - T™, (A20) 



satisfies the equation 



which explicitly links the antisymmetric part of the energy-momentum tensor (jA7J) with the antisymmetric tensor 
w a[y<r]_ when energy-momentum is conserved (i.e., d^T^ ee for v = 0, 3), the total angular-momentum (|A16[) is 
also conserved, while the orbital and spin angular momenta satisfy the equation 
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Hence, the angular-momentum conservation law d^J^ 01 ^ = relates the orbital angular-momentum density L°^ a ^ 
with the spin angular- momentum density S°^ a ^ in a form that is similar to Eq. (| 1 19[) . 

Lastly, we note that the canonical energy-momentum tensor T^ v can also be symmetrized (34|-[36j directly by adding 
the term d a Y\ a ^ 1 ' (where Y^ a ^ 1 ' = — £^ IQ ] l/ ) and require that the new energy-momentum tensor be symmetric: 

ee T^ V + d a T,^ v ee T Vfi , (A22) 

while satisfying the same energy-momentum conservation law d^T ee d^T^ (since df ia ^ a ^ l/ ee 0). The Belinfante- 
Rosenfeld solution 37] to Eq. (|A20I) is expressed as 

ji[an]u = }_ ^gu[an] _|_ gfJ,[av] _ ga[fj,u]^ (A23) 

so that Y\ a ^ 1 ' — El 01 "]* 4 = — S a ^ v \ Hence, the spin angular-momentum density S a ^ v \ which is driven by the 
asymmetry of the canonical energy-momentum tensor [Eq. (|A20p ]. can also be used to obtain a symmetric energy- 
momentum tensor dA"22]) . where T^ a ^ v is defined by Eq. (|X23]> . 
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